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Abstract 

At zero temperature, the 3-state antiferromagnetic Potts model on a square lattice 
maps exactly onto a point of the 6-vertex model whose long-distance behavior is 
equivalent to that of a free scalar boson. We point out that at nonzero temper- 
ature there are two distinct types of excitation: vortices, which are relevant with 
renormalization-group eigenvalue ^; and non- vortex unsatisfied bonds, which are 
strictly marginal and serve only to renormalize the stiffness coefficient of the un- 
derlying free boson. Together these excitations lead to an unusual form for the cor- 
rections to scaling: for example, the correlation length diverges as = J/kT — > oo 
according to £ ~ Ae 2 @(l + bf3e~@ + • • •), where b is a nonuniversal constant that 
may nevertheless be determined independently. A similar result holds for the stag- 
gered susceptibility. These results are shown to be consistent with the anomalous 
behavior found in the Monte Carlo simulations of Ferreira and Sokal. 



1 Introduction 



Corrections to scaling at continuous phase transitions have long been understood 
within the framework of the renormalization group [HJ. Irrelevant operators give rise 
to power-law corrections to scaling in, for example, the temperature-dependence of the 
correlation length, of the form 

£ ~ A\t\-"(1 + C\t\» + ...) (1.1) 

where u > is a universal correction-to-scaling exponent, and A and C are nonuniversal 
amplitudes. (For a critical point at temperature T c ^ 0, the scaling variable t is con- 
ventionally defined by t oc T — T c .) In the case of a marginally irrelevant operator, both 
multiplicative and additive logarithmic corrections can occur, typically of the form 

£ ~ A^ilog^- 1 ) 

However, it is possible for more exotic dependences to arise. This paper considers just 
such an example, namely the antiferromagnetic 3-state Potts model on a square lattice. 
This model has the Hamiltonian 

(ij) 

where the sum is over nearest-neighbor pairs of vertices on the square lattice, at each 
of which is a degree of freedom s« taking one of three possible states. The coupling is 
antiferromagnetic, that is, J > 0. 

One reason that this model is exotic is that its critical point occurs at T = 0, |3|, [5], 
H) S H) H- Nevertheless, this is a bona fide critical point, exhibiting, for example, power- 
law decay of the order-parameter correlation functions. At nonzero temperature, these 
correlation functions decay exponentially with a finite correlation length £. In the past 
there has been some confusion concerning the correct choice of scaling variable at such 
zero-temperature critical points; it does not necessarily correspond to setting T c to zero 
in the above definition of t (i.e. to taking t oc T). In fact, this question is answered quite 
explicitly, at least for non-quantum critical points, by the renormalization group: the space 
of scaling variables is the tangent space at the fixed point to the manifold parametrized 
by the Boltzmann weights of the model in question. In our case, the low-temperature 
configurations correspond to modifying one of the ground states by allowing a nonzero 
density of nearest-neighbor bonds where s, = Sj (we henceforth call these unsatisfied 
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bonds). Each unsatisfied bond occurs with Boltzmann weight e _/3 (where (3 = J/kT), and 
therefore the correct scaling variable t should be linear in e _/3 . 

The other, more interesting, reason that this model is unusual is that, as we shall dis- 
cuss in detail in the next section, switching on the temperature excites not just one, but 
two, scaling operators. One of these operators turns out to be relevant, with renormalization- 
group eigenvalue y = | (hence scaling dimension x = 2 — y = |), corresponding to 
v = 1/y = 2, so that the leading behavior of the correlation length is £ oc (e~P)~ u = e 2f3 . 
The other operator is marginal — not marginally irrelevant, but strictly marginal in the 
sense that, taken alone, it would generate a line of fixed points with continuously varying 
exponents. The main result of this marginal operator is to give the exponent v in the 
foregoing expression an effective dependence on the scaling variable w = e _/3 itself: thus 
we find that 

f ~ A{e~ p )- V ^ B) ~ Ae 2 \\ + b/3e~ p + • ■ •) (1.4) 

where b = dv e g(w) / dw\ w= o is presumably nonuniversal. It turns out, however, that in this 
model we can estimate the value of b by an independent calculation. In particular, simple 
qualitative arguments show that b < 0, so that the asymptotic value of £,/e 2/3 should be 
reached from below as (3 — > oo (and likewise for the staggered susceptibility). 

Such an increase, following a minimum, was in fact found in extensive Monte Carlo 
simulations of this model by Ferreira and Sokal [TIJ. These authors attempted to fit their 
data with a variety of "standard" forms of corrections to scaling, including logarithms. 
Though some of these Ansatze gave reasonable fits, none had a plausible theoretical basis 
1T0| , Section 7.1] for detailed discussion). As we shall show in Section |, the anomalous 
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behavior found by Ferreira and Sokal is consistent with the unusual form ( |L4| ) over its 
expected range of validity, when b is given the value that we extract independently in 
Section || 

Our ability to make such a theoretical analysis depends on a mapping of the model 
at zero temperature to a discrete height model that can be connected (via a nonrigorous 
but convincing renormalizat ion-group argument) to the continuum theory of a free scalar 



boson [g, [11], |9fl. It is within the height-model approach that we are able to disentangle 
the behavior at nonzero temperature and argue that the two types of excitation, leading 
to the unusual scaling form (El), are present. 

The layout of this paper is as follows: In Section ^| we discuss the mapping to a height 
model and make the important observation that in this model (unlike some other models 
with zero-temperature critical points) the height mapping is well-defined also at nonzero 
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temperatures. It is then straightforward to identify the two types of excitation in height- 
model language. The most relevant excitation has the nature of a vortex within the free- 
boson description, and standard Coulomb-gas methods then lead to the prediction that 
v — 2. The other type of excitation has zero vorticity and, we argue, merely renormalizes 
the compactification radius (or stiffness coefficient) of the free field. However, this is not a 
rigorous argument, and it is important to check it independently. This we do in Section [3], 
by a direct numerical investigation of the modified height model in which only excitations 
of the second type are allowed, that is, the vortices are suppressed. In this case, the height 
model may be represented as a modified vertex model, which we analyze using transfer 
matrices and finite-size scaling. Our results show clearly that these excitations are indeed 
strictly marginal: the modified vertex model is still critical and continues to have the 
central charge c = 1 characteristic of a free boson. We directly measure the dependence 
of the stiffness coefficient on the fugacity w = e _/3 of the excitations. In Section f| we put 
these pieces of information together to predict how the effective critical exponents like 
v e g depend on e~^, and hence extract the value of the parameter b in Eq. ( |1.4j ) for both 
the correlation length and the susceptibility. Finally we compare these predictions with 
the Monte Carlo data of Ferreira and Sokal [|10|| . In the Appendix we show how various 
non-universal quantities arising from the first-order effects of the non- vortex defects may 
be related to one another. 



2 Height model and low-temperature excitations 

It is convenient to label the Potts states Sj by the integers 0, 1, 2. The vertices i of 
the square lattice are labeled by integer- valued coordinates (mj,ni), and the lattice is 
divided into even and odd sublattices on which rrii + n« is even or odd, respectively. Let 
us introduce the variable 77, = ~[1 — (— l) m »+ ri i] ! which takes the values or 1 according to 
whether i is on the even or the odd sublattice. We then define height variables hi G Z/6Z 
by 

hi = 2 Si + 3r)i (mod 6) , (2.1) 
as illustrated in the table below: 



Potts state 


Even Subl. 


Odd Subl. 
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This gives a 1-1 correspondence between configurations of the 3-state Potts model and 
configurations of the height model satisfying the constraint that hi is even (resp. odd) 
whenever i is on the even (resp. odd) sublattice. This latter constraint may equivalently 
be imposed by fixing the height at the origin to be even and further demanding that 
heights on neighboring vertices i and j satisfy 

\hi -hj\ = 1 or 3 (mod 6) . (2.2) 

We then have S, ^ Sj if and only if — hj\ — 1 (mod 6), and S, = Sj if and only if 
\hi — hj\ =3 (mod 6). 

Let us consider first the model at zero temperature. Then the only allowed configu- 
rations are antiferromagnetic ground states (sj 7^ Sj for all nearest-neighbor pairs i, j); in 
the height model this corresponds to replacing ( [2.2|) by the more restrictive condition 

\hi-hj\ = 1 (mod 6) . (2.3) 

It then follows that the height field hi G Z/6Z can be "lifted" to a height field hi G Z 
so that now the height difference across an edge is ±1 tout court (not just mod 6): self- 
consistency is ensured (at least with free boundary conditions) by noting that the change 
Ah around any plaquette will be zero (if four numbers ±1 add up to mod 6, they 
must necessarily be two +l's and two — l's, hence add up to 0). Furthermore, this lifting 
is unique up to an overall shift by a multiple of 6. In cylindrical or toroidal boundary 
conditions, the height field hi might fail to be globally well-defined (i.e. it might have a 
nonzero tilt), but the gradient field V/i is still well-defined and curl-free. 

The restricted height model (|2.3| ) can be equivalently be mapped onto an arrow model 
on the dual lattice: we assign to each edge e of the dual lattice, dual to the edge ij of 
the original lattice, the orientation obtained by a +90° turn from the direction in which 
the height change is +1 (mod 6). This arrow field is simply the dual of the vector field 
V/i; it is therefore conserved at each vertex of the dual lattice. Each vertex of the dual 
lattice thus has precisely two inward-pointing arrows and two outward-pointing arrows, 
so that the allowed configurations of the arrows at each vertex are precisely those of the 



6- vertex model [13]. Moreover, the Boltzmann weights are those of the symmetric point, 



i.e. all weights equal. It is well known from the exact solution of the 6-vertex model 
that this point is critical. 



As was pointed out by Henley [0, 11 , this criticality is quite easily understood from 



the point of view of the renormalization group applied to the lifted height model.Q It 
1 One might ask why the RG has to be applied to the lifted (Z-valued) height model and not to 
the original (Z/6Z-valued) height model. The answer is that the coarse-graining process necessarily 
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is reasonable to guess that the long-wavelength behavior of the lifted height model is 
controlled by an effective coarse-grained Hamiltonian of the form 

K 



Sa = 



^-{dhf - Acos(2tt/i) 



(2.4) 



where K is the stiffness constant. The gradient term in (|2.4j ) takes into account the 



entropy of small fluctuations around Henley's "ideal states" |L2|, [L3|, §, [□], [|; the second 
term is the so-called locking potential, which favors the heights to take their values in 
Z. We then expect that there exists some constant K T such that for K < K r (resp. 
K > K T ) the locking potential is irrelevant (resp. relevant) in the renormalization-group 
sense. Thus, if K < K r our surface model is "rough" and its long-wavelength behavior 
can be described by a massless Gaussian model: 

[h(x) - h(y)] 2 ) ~ -i-loglar-yl (2.5) 



7lK 

for \x — y\ ^> 1; in this case, the original zero-temperature spin system is critical, and all 
its critical exponents can be determined in terms of the single constant K. In particular, 
the scaling dimensions of spin-wave (or vertex) operators e iah ^ are given by 

x a = a 2 /4nK. (2.6) 



From (|2.1| ) the staggered order parameter (_i) m *+ n ; e 27rs i/ 3 corresponds to a = tt/3, so that 



its correlation function decays with an exponent //stagg = 2%tt/3 = tt/18K. By comparison 
with the exact result 77 stag g = | from the 6- vertex model 0, 0, it follows that K must 
take the value 7r/6. At this value, the locking potential cos(27r/i) has scaling dimension 
x 2lT = 6, so that it is indeed highly irrelevant. By the usual scaling law we then obtain 
the susceptibility exponent (7/V) s ta g g = 2 — r/stagg = 5/3.0 

Let us now consider the model at nonzero temperature: note that the height model 



( |2.1| ) / ( OP continues to give a complete description. Nonzero temperature amounts to 
allowing the height difference between neighboring sites to take the value 3 (mod 6) as 
well as ±1 (mod 6): the former correspond to unsatisfied bonds. Each unsatisfied bond 
will be weighted by e _/3 , so that, at low temperature, the unsatisfied bonds will be very 
dilute. Let us therefore consider first the effect of just one isolated unsatisfied bond. 
This bond belongs to two neighboring plaquettes, and there are two cases to consider, 



introduces fractional weights, which would cause ambiguity in R/6Z but which give rise to a well-defined 
averaging operator in R. 

2 This value has been confirmed numerically by several authors (T^, 16, |[ [To| . 
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according as the vorticities on those two plaquettes have (a) opposite signs or (b) the 
same sign: 

2 1 

s = — s 
2 1 

(a) 

where the unsatisfied bond is indicated with a dash 
rations mod 6 are 

1 2 

h = — 3 h 
1 2 

(a) 

(here we have taken the center-left vertex to lie on the even sublattice). Let us now ask 
whether these heights hi 6 Z/6Z can be lifted to heights hi 6 Z while satisfying the 
condition 

\hi-hj\ = 1 or 3 . (2.9) 

This amounts to asking whether the bond in ( |2.8| ) with \hi — hj\ =3 can be assigned a 
sign consistent with that of the two adjacent plaquettes. In case (a), the answer is yes; in 
case (b), the answer is no. We therefore call situation (a) a non-vortex unsatisfied bond, 
and situation (b) a vortex (of strength 6). 

This situation can alternatively be viewed in the 6- vertex picture (Figure []]). The 
arrow dual to the unsatisfied bond may be thought of as an arrow of strength 3, but 
of indeterminate sign. There are then two possible ways in which this defect may be 
healed locally: either we can choose a definite orientation for the triple arrow, and impose 
strict conservation at the neighboring vertices, corresponding to the situation shown in 
Figure |l]a; or we can relax this condition and allow a violation of arrow conservation at 
one of the vertices, as in Figure |l]b. In the latter case there is a net arrow flux of ±6 
out of the region of the defect, and this will persist to larger distances if there is strict 
conservation elsewhere. 

Clearly this kind of defect has a topological nature: it corresponds to a vortex (or 
antivortex) in which the height field changes by ±6 in encircling the defect. Outside the 
core of the defect, the continuum description ( |2.4| ) should still apply. Thus we may write 
h(r) « (±6/2tc)6 + h f , where 9 is the polar angle and h' has zero vorticity. Substituting 
into (|2.4j ) we then find that the defect has additional reduced free energy (i.e. negative 



2 1 

= — 

(2.7) 

12 v ' 

(b) 

The corresponding height configu- 

1 2 

= — 3 

(2.8) 

5 4 v ; 

(b) 
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entropy) 

\ K {^) J^/r 2 ~ (9K/n)\og(R/a) (2.10) 

where we have introduced a short-distance cutoff a, and R is the size of the system. Thus 
a single isolated vortex has zero probability of occurring. However, configurations contain- 
ing several vortices and antivortices with zero total vorticity have a nonzero probability 



as R — > oo. The scaling dimension xy of a vortex operator may be determined §L7[ p. 121] 
from the power-law decay of the vortex-antivortex correlation function (i.e. the partition 
function with a vortex-antivortex pair introduced) or, more easily, directly from ( 2.10|) as 



x v = 9K/n = | (2.11) 

at K = 7r/6. If the vortices are the most relevant perturbation at nonzero temperature, 
which we shall argue is the case, then (|2.11| ) leads to the prediction v = 1/yy = 1/(2 — 
Xy) = 2. Since the scaling variable t is proportional to the fugacity e _/3 of the vortex, we 
predict that the correlation length diverges as £ oc {er^)~ v = e 2/3 as ft — > oo. 

On the other hand, the non-vortex defect shown in Figure [l|a has finite additional 
negative entropy. It may be viewed as a tightly bound pair of a vortex and antivortex of 
strength ±3. Thus, at large distances, it corresponds to a dipole height field 

h{r)~D?-^ (2.12) 

where n is a unit lattice vector, and D is a non-universal constant whose value is difficult 
to determine analytically, since the size of the vortex-antivortex pair is of the order of the 
lattice spacing, at which scale the continuum action ( |2.4| ) is inapplicable. (Nevertheless, 
in the Appendix we show how the value of D may be related to other non-universal 
numbers that we have measured directly.) The dipole field ( f2.12j ) gives an infrared-finite 
but nonzero contribution to Sq (of course the integral should still be cut off at short 
distances: it diverges as a~ 2 and so is strongly dependent on the precise form of the 
cutoff). The numerical value of this nonuniversal constant will be measured in Section [3[ 
As already noted, the non-vortex defect has the property — at least in the simple 
case when it is isolated — that it is possible to assign a definite value +3 or —3 to the 
height change Ah along the defect edge in such a way that the gradient field Vh remains 
curl-free, i.e. such that the height variables are lifted locally to Z. (In the 6- vertex picture, 
this means that it is possible to assign a definite arrow value +3 or —3 to the defect edge 
in such a way that strict conservation holds.) Even when the defects are not isolated, 
but form finite clusters, as long as there is no net flux of arrows from each cluster, it 
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should still be possible to lift the height variables locally to Z in a unique manner, except 
possibly in the core of the cluster. (An example of ambiguity within the core of the cluster 
is shown in Figure 0.) We therefore expect that non- vortex defects do not disturb the 
renormalization-group flow of the model towards a free boson theory as in fl2.4|) . On the 
other hand, there is no reason to suppose that a model with a finite density of defects will 
renormalize onto the same value of the stiffness constant K as is obtained in the absence 
of defects; rather, we expect that K will depend on the fugacity w = of the non- vortex 
defects. It is the purpose of the next section to test this hypothesis numerically and to 
estimate the dependence K e g(w) in the model with vortices suppressed. 

In general, the renormalization-group flows in a model with a scale-dependent stiffness 
constant K(£) and vortex fugacity y(£) are well-known |]18| to be (for a vortex of strength 

6) 

dy/d£ = (2-9K/Tr)y + 0(y 3 ) (2.13a) 
dK^/di = Ay 2 + 0{y 4 ) (2.13b) 

The coefficient of y in the first equation reflects the fact that our vortices have scaling 
dimension xy = 9K/ir. The coefficient A in the second equation is nonuniversal (since 
it depends on the normalization of y) but of order unity; it is independent of K because 
the vortex-vortex interaction is oc K, and the screening term from tightly bound pairs is 
oc K 2 , so that dK/d£ oc K 2 y 2 . 

Defining u = tt/6 — K, the RG flow (2.13) may be rewritten as 

dy/d£ = (| + 9u/Tt)y + 0(y 3 ) (2.14a) 
du/d£ = A'y 2 + 0{y\y 2 u) (2.14b) 

where A' = (tt/Q) 2 A. Away from the point where u = — 7r/18 (where the vortices become 
marginal), these equations may be rewritten in the standard way [19| in terms of non- 
linear scaling variables u = u + 0(y 2 ) and y = y + 0(y 3 ) so that the flow equations take 
the simple form 

dy/d£ = (± + 9u/n)y (2.15a) 
du/d£ = (2.15b) 

with no higher-order terms. These equations are to be integrated with the initial condi- 
tions y(0) =w = e~ 13 and K(0) = K cS (w), hence y(0) = y(0) + O(y(0) 3 ) = w + 0(w 3 ) and 
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fi(0) = u(0) + O(y(0) 2 ) = f - A^H + 0(w 2 ) = 
and 

y(^) = [ w + 0(w 3 )} exp [(i + 9m(0)/tt)£ 



AT^(0)w + 0(w 2 ). Thus w(£) = fi(O) 



(2.16) 



The correlation length satisfies the homogeneous renormalization-group equation 



(2.17) 



Assuming in the standard way that £ = £ = 0(1) when y{£) = 0(1) [i.e., when y(£) 
0(1)] then gives the prediction 



~ £o [w + O(^ 3 )]- [2+(36/w) ^« (0)lu+o(? " 2)] 

= £ exp [2/5 + &/3e _/3 + 0(/3e" 2/3 )" 

= £ e 2/3 [l + bpe-P + \b 2 l3 2 e~ 213 + 0(/3 e - 2/3 ) 



(2.18a) 
(2.18b) 
(2.18c) 
(2.18d) 



where b = (36/7r)i^ ff (0). In writing the above we have been careful to show where the 
neglected higher-order terms enter. Simple qualitative arguments (see Section |j) show 
that K' eS (0) < 0, so this implies that the asymptotic value of £(/3)/e 2/3 is attained from 
below. 

Next we turn to the staggered susceptibility. A staggered field h stagg satisfies the 
renormalization-group equation 



dh stagg /d£ = (2 - x 7r/3 )/i stagg + 0(h 3 st ) 



(2.19) 



where, as discussed above, = n/36K. As usual, the singular part of the reduced free 
energy per unit area transforms according to f(u(0),y(0), h stagg (0)) = e~ 2e f(u(£),y(£), h staug , 
so that the staggered susceptibility Xstagg ~ d 2 f / dh 2 tSugg satisfies 



Xstagg(«(0),y(0)) = exp 



where the integral in the exponential is 



[2-2x w/3 (£')]d£' 



Xstagg(w(^) S y(^)) 



^2.20) 



[| - (2w/tt) + 0(y(£' ) 2 )} d£' ~ (§ - 2u/tx)£ + const + 0(w 2 ) . (2.21) 



Choosing £ as before and assuming that at this scale Xstagg = Xo — 0(1) then gives 

r 1 r\t 2\i-(|-2u/7r)/(i+9?i/7r)+0(«) 2 ) 



Xsta; 



gg 



(2.22) 
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where, once again, u = —K' cS (0)w + 0{w 2 ). Equation (|2.22| ) then simplifies to the final 
result 



Xst agg (/3) ~ Xoex P \f(3 + b'f3e^ + 0(f3e-^) 



1 + tipe-P + \b' 2 p 2 e~ w + 0(/3e 



(2.23a) 
(2.23b) 



where b' = (64/ir)K' e{[ (0) = ^-b. Thus the asymptotic limit of Xstagg(/3)/e 3 13 is also 
attained from below. 



3 Extracting the stiffness constant 

In this section we analyze a model in which the vortex excitations present in the full 
model are deliberately suppressed, in order to understand the effects of the (less relevant) 
non-vortex defects like those shown in Figure |l]a. As discussed in the previous section, 
this may be achieved in the height model by allowing height changes Ah of ±3 as well 
as ±1, while continuing to insist that the gradient field V/i be strictly curl- free. Edges 
with Ah = ±3 will be assigned a fugacity w = e" 13 . Equivalently, we may work in the 
dual arrow model: edges in this model are allowed to carry flux ±3 as well as ±1, and we 
impose strict flux conservation at the vertices. The partition function of the arrow model 
is therefore 

Z = 5>"», (3.1) 
Q 

where N 3 is the number of triple arrows, and the sum J2g extends over all possible 
flux-conserving configurations of the 44-vertex model just defined.^ We wish to test the 
hypothesis that the continuum limit of the statistical model defined by ( |3.1|) is a free-field 
theory with action given by (|2.4j) , with a stiffness constant K that depends on w. In this 
section we shall see how it is possible to verify this hypothesis numerically, and at the 



3 Barbero et al. [£(J recently studied a closely related 14-vertex model, in which zero or one edges 
of strength ±3 (but not two or four such edges) are allowed to be incident on each vertex. For small 
fugacity w, this model should be essentially equivalent to our model; in particular, we expect it to remain 



in the "rough" phase for small w, with a continuum limit given by the massless free field (2.4) with a 
stiffness constant K that depends on w. However, this model cannot stay critical for all w > 0, since 
in the limit w — > oo it freezes into one of a finite number of ground states (for example, vertices 11 and 
16 of ||^, Figure 1] on even and odd sublattices, respectively), ft would be interesting to understand 
what happens in-between. Unfortunately, the study of Barbero et al. is limited to the region < w < 1 
(corresponding to fie^ > in their notation; we are at the symmetric point /3e = 0). We thank Giorgio 
Mazzeo for bringing this paper to our attention. 
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same time extract quite accurate values of K e g(w) as well as its first few derivatives at 
w = 0. 

Imagine defining the model ( |3.1| ) on a cylinder of circumference L, with periodic 
boundary conditions in the transverse direction. Clearly, due to the flux conservation it 
splits up in a direct sum of theories with a fixed net arrow current Q in the longitudinal 
direction. The corresponding height field h(i) = h(x, t) is multiple- valued in the transverse 
coordinate x, with "tilt" Q: h(x + L,t) = h(x,t) + Q. Such a constant tilt can, however, 
easily be gauged away by setting 

h( x ,t) = Qx/L + h(x,t) , (3.2) 

where now h(x + L,t) = h(x,t). For Q/L ^ 1 the field h can be assumed to describe 
the Q = sector of the original model. Therefore, inserting ( |3.2| ) into ( |2.4| ), we see that 
a nonzero current Q simply shifts the action (free energy) per unit area by an amount 

A/(Q) = ^Mg. (3.3) 

This is in turn related to the shift in the logarithm of the largest eigenvalue of the corre- 
sponding transfer matrix, which can then be used to check the above hypothesis and to 
extract numerical values of K eS (w). 

To this end we have constructed the transfer matrix of the model ( |3.1|) in the 4 L - 
dimensional basis associated with all possible arrow configurations for a layer of L vertical 
bonds. As mentioned above, the transfer matrix is block-diagonal according to the value 
of Q. Furthermore, Q must have the same parity as L. Using standard sparse-matrix 
techniques we have been able to diagonalize the various sectors of this matrix for widths 
up to L max = 10. Due to parity effects we limit the discussion to even L in the following, 
although odd L yield compatible results. 

To check that our model is indeed described by a Gaussian theory for all w > 0, we 
begin by examining the central charge. This can be extracted from the finite-size scaling 



of the free energy per unit area [21 



/o(i) = /oN-^ + '- , (3-4) 

where fo(L) = — logAo(-fz) and Ao(£) is the largest eigenvalue of the transfer matrix in 
the Q = sector. As usual, the convergence can be sped up by including a nonuniversal 
1 /L 4 correction, so that three consecutive system sizes are needed to fit the above formula 
. The resulting finite-size estimates for c as a function of w are shown in Table [I]. As 
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w 


c(6) 


c(8) 


c(10) 


0.0 


0.9751 


0.9885 


0.9958 


0.5 


1.0079 


0.9971 


0.9990 


1.0 


0.9430 


0.9812 


0.9968 


1.5 


1.3336 


1.0448 


0.9407 


2.0 


1.6105 


1.3984 


1.2046 


3.0 


0.9486 


0.9335 


0.9771 


10.0 


0.9618 


0.9882 





Table 1: Effective central charge c as a function of the Boltzmann weight w. The three- 
point fits to ( |3.4| ) based on system sizes L, L — 2 and L — 4 are labeled as c(L). 

expected, they strongly suggest the ^-independent value c = 1 in the L — > oo limit, thus 
corroborating the analytical arguments given in Section |2]. (For reasons that we do not 
fully understand, the finite-size effects are particularly strong at w ~ 2. This would be 
an interesting question to investigate further.) 

We now turn our attention to the extraction of the stiffness constant K e q(w). For 
each L < L max we have verified that the Q-dependence of the free energy is indeed of 
the form ( |3.3|) for Q ^ L. However, as the quadratic behavior must (and does) saturate 
when Q becomes comparable to L, we have based our estimates of K e g (w) on only the two 
lowest permissible values of Q. For even L these values must be even, due to the parity 
observation made above; this suggests that we focus on Q = and Q = 2. But at this 
point special attention must be paid to the limit w — > oo, in which only triple arrows are 
allowed, so that the model reduces to the usual 6-vertex model with triple-size heights. 
Such configurations are incompatible with a total flux of Q = 2. Thus, for generic (i.e. 
non-small) values of w we should restrict attention to values of Q that are multiples of 6, 
e.g. Q = and Q = 6. 

The values of K e g(w) extracted by fitting ( |3.3|) at Q = and Q = 6 are displayed in 
Table |2] for selected values of w and L. We observe a leading finite-size correction in 1/L 2 , 
reflecting the 1/L 4 correction present in c. Thus, fitting each pair of successive system 
sizes to the form K e s(L) = K e g(oo) — const/L 2 , we obtain two independent estimates 
of i^eff(oo); and judging the error bar from the small residual size dependence of these 
estimates we obtain the final result shown in the rightmost column of Table 0. For w = 
we find excellent agreement with the exact result for the 6-vertex model, K e g(0) = 7r/6. 
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w 


K(L = 6) 


K(L = 8) 


K(L = 10) 


Extrapolation 


0.0 


0.505099 


0.512823 


0.516574 


0.52359(1) 


« tt/6.000 


0.5 


0.207033 


0.205920 


0.205408 


0.20450(1) 


« tt/15.4 


1.0 


0.125262 


0.125426 


0.125496 


0.12562(2) 


~ tt/25.0 


1.5 


0.091195 


0.092921 


0.093571 


0.0947(4) p 


d tt/33.2 


2.0 


0.064168 


0.064762 


0.064972 


0.0654(2) s 


n tt/48.1 


3.0 


0.057232 


0.058096 


0.058511 


0.05925(5) 


« tt/53.0 



Table 2: Finite-size estimates for the stiffness constant K e g(w) as well as their extrapo- 
lation to L = oo. 



Quantity 


L = 4 


L = 6 


L = 8 


L = 10 


Extrapolation 


K'(0) 


-0.330292 


-0.452044 


-0.502110 


-0.527432 


-0.58 (1) 


K"(0) 


-3.10237 


-2.69860 


-2.49302 


-2.37655 


-2.1 (1) 



Table 3: Finite-size estimates and extrapolations for the derivatives dK/dw\ w=0 and 
d 2 K/dw 2 



\w=0- 



For w — > oo the triple-height model should have correlations ( |2.5|) nine times as large, 
hence a stiffness constant one-ninth as large, i.e. K e g(oo) = 7r/54; this too is confirmed 
numerically. In-between, K e s(w) is a monotonically decreasing function of w; this makes 
sense heuristically, because allowing Ah = ±3 increases the variance of the distribution of 
nearest-neighbor height differences, which should lead (or so one naturally expects) to a 
larger variance also for the long-distance height differences, which by ( |2.5|) are proportional 
to K-\ 

Exactly at w = the above remark that the Q = 2 sector introduces frustration does 
not apply. The evaluation of K and its first few derivatives at w = can therefore be 
based on the Q = and Q = 2 sectors, yielding a better precision. To perform the 



derivatives, we simply numerically differentiate the free energies ( |3.3|) . The results for 
K'(0) and K"(0) are given in Table §. 

Our transfer-matrix results can also be used to compute the average energy density, 
given here simply as the probability of having a triple arrow, e(w) = (N 3 )/N. Using (|3.1| ) 



we have e(w) = —wdf/dw, and based on the first two derivatives of the free energy we 
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find 

e(w) = 0.2179(2)w + 0.696(l)w 2 + O(«; 3 ) (3.5) 

in the model with vortices suppressed. This can be compared with the Monte Carlo 
results of Ferreira and Sokal for the full 3-state Potts antiferromagnet (i.e. with vortices 
included), based on a fit to the low-temperature data 0, Section 4.3 and Figure 13]: 

epottsH = 0.21777w + 1.65303w 2 + O(u; 3 ) . (3.6) 

The fact that the leading terms of the two expansions agree could easily have been antic- 
ipated, since vortices necessarily come in pairs, and thus enter only at order 0{w 2 ). 

4 Comparison with Monte Carlo results 

Putting together the results of Sections ^ and |, we deduce the theoretical predictions 

log £03) = 2(3 + log£o + b(3e- p + 0(pe~ 213 ) (4.1) 

\og X {P) = f P + logxo + b'Pe-P + 0(/3e" 2/3 ) (4.2) 

where the nonuniversal constants b and b' take the values 

36 

b = —K' eS (0) w -6.65(11) (4.3) 

7T 

b' = -K s (0) « -11.82(20) (4.4) 

71 

and we write \ as a shorthand for Xstagg- We are now ready to compare these predictions 
with the Monte Carlo data of Ferreira and Sokal [JTU], Table 4, L min = 128]. 

Let us begin with the correlation length. In Figure ||] we plot £ (/3) / e 2/3 versus /3; the 
rise at ft ft 3.4 (corresponding to £ <; 75) is seen clearly. In Figure f| we replot the same 
data as log£(/3) — 2/3 versus Pe~^\ there is a lot of curvature, and it is unfeasible to extract 
reliable estimates of the limiting slope, but the data for /?e _/3 < 0.031 (corresponding to 
P <; 5.1 and £ ^ 2600) are at least compatible with the predicted slope —6.65, provided 
that we choose an intercept log £o = —2.11(2). In Figure||we plot log £(/?) — 2/5 + 6. 65/5e _/3 
versus Pe~^\ for modest values of P the corrections to scaling are actually stronger than 
in Figure f| (!), but for large P they are weaker and it is at least plausible that the 
curve is asymptotically horizontal. Finally, in Figure |6] we plot log£(/3) — 2/? + 6.65/?e _/3 
versus /3e~ 2/3 ; now the curve is compatible with linearity over the much wider range 
Pe~ 213 < 0.0019 (corresponding to P > 3.8 or £ > 170), and we estimate an intercept 
l°g£o — —2.12(1) and an asymptotic slope C = 123(4). 
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Next we analyze the staggered susceptibility. We omit the plot of x(P) I e^ 10 / 3 - 1 ^ versus 
(3; it looks a lot like Figure |3|. In Figure |7| we plot logx(P) — y/3 versus /3e _/3 ; once 
again there is a lot of curvature, and it is unfeasible to extract reliable estimates of the 
limiting slope, but the data for /3e _/3 < 0.031 are compatible with the predicted slope 
—11.82, provided that we choose an intercept logxo = —2.52(3). In Figure || we plot 
logx(P) — y/3 + 11.82/3e _/3 versus /3e _/3 ; for modest values of (3 the corrections to scaling 
are once again stronger than in Figure |7|, but for large /3 they are weaker and it is at 
least plausible that the curve is asymptotically horizontal. Finally, in Figure |] we plot 
logx(P) — y@ + H-82/5e _/3 versus /3e _2/3 ; now the curve is compatible with linearity over 
the much wider range (3e~ 2lS ^ 0.0016 (corresponding to (3 ^ 3.9 or £ > 210), and we 
estimate an intercept logxo — —2.54(1) and an asymptotic slope C = 232(8). 

In conclusion, the available Monte Carlo data are compatible with the theoretical 
predictions (|4.lQ / (|4.2| ), although the evidence for these predictions — and in particular 



for the predicted values ( |4.3|) / (4~4) of the nonuniversal constants b and b' — is admittedly 
less than overwhelming. Our analysis does, in any case, give a simple explanation of why 
the limiting values of £(/?)/e 2/3 and x(P) I 'e^ 10 / 3 -^ are approached from below. It would 
be useful to obtain higher-precision Monte Carlo data at correlation lengths £ ^ 1000 
(corresponding to (3 ^ 4.6) in order to make a better test of our theoretical predictions. 
It is curious that the asymptotic behavior in this model is attained only at rather large 
correlation lengths. 
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A Relations between non-universal constants 



In this appendix we show how dK/dw\ w=0 is related to two other non-universal con- 
stants: D in (gig ) and d = de/dw|^ =0 in (TO /(TO- 

In the Coulomb gas picture of vortices interacting with a strength oc K, the tightly- 
bound vortex-antivortex pairs of vorticity ±3 (Figure [l|a) act to screen this interaction, 
thus decreasing the effective K. A very similar effect happens in the Coulomb gas picture 
of the XY model, and we show that a calculation similar to that employed in the Kosterlitz 
renormalization-group approach may be used here to estimate the shift in K, to first 



order in w = e~@. 

Rather than evaluating the interaction between two vortices, it is simpler to consider 
the correction to the height correlation function, which, in the absence of defects, has the 
form 

([/i(n) - h(r 2 )] 2 ) ~ (TrlO -1 log 1^-^1/0 + const. (A- 1 ) 
Introducing defects of the form 



W) = ^ V""^ (A.2) 
r - r 2 



where n is a unit lattice vector, the first-order correction to the (h(ri)h(r 2 )) term in (|A.1| ) 



is 

E _D 2 £ (n "| ri " r0 ? 2 ) | n " (r2 7 2 r ° )) (A.3) 



e -P e -s M 



ro,n 



ri - r r r 2 - r 



where S^et is the action (relative negative entropy) of the defect. Summing over the four 
orientations of n turns the numerator into 2(r% — r ) • (r 2 — r ). In the continuum limit, 
( |A.3|) may be written 

2e-Pe- s *«D* [ ^ ~ T °] ' ^ - r ° } dh (A.4) 
J \ti - r \ 2 |r 2 - r | 2 

Apart from pieces which contribute to terms independent of |ri — r 2 |, the numerator in 
( |A.4| ) may be rewritten as — |[(ri — r ) — (r 2 — r )] 2 = — |(ri — r 2 ) 2 , so that ( |A.4| ) becomes 



_ e -/3 e -s dC f D 2 f l r i r 2l d 2 ro e -^ e ~ 5dcf D 2 -47rlog|ri -r 2 |/a (A.5) 

J ri-r M r 2 -r r 



(This last result may be most easily seen by evaluating the dependence on the cut-off a 
close to ro ~ ri and ro ~ r 2 .) From (|A.5|) may be read off the first-order correction to 
the stiffness constant 

5K/K 2 = -2tt ■ D 2 ■ 4 7 re- /3 e- 5dcf . (A.6) 
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Setting K = n/6 thus gives 

8K = -(2tt 4 /9) • D 2 ■ e- p e~ s ^ . (A.7) 

On the other hand, the first-order correction of such defects to the free energy per 
vertex is 

Sf = -4e- /3 e" Sdcf = -cie"" (A.8) 



where c\ is the number appearing in (3J5) / (|3.6| ) and numerically determined to be ~ 0.218. 
This gives the value dK/dw\ w= o = — (7r 4 /18)D 2 ci, leading to the estimate D pa 0.70. This 
is in rough agreement with various rather crude estimates which may be made on the 
assumption that the effective Hamiltonian Sq is valid down to the lattice scale. 
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Figure 1: Examples of the two types of defect: (a) non- vortex unsatisfied bond, and (b) 
vortex. In both cases, two neighboring sites (labeled here by heights and 3) are in the 
same Potts state (here s = 0). Case (a) is a local defect with finite relative entropy, while 
(b) corresponds to a vortex in the 6- vertex model, with logarithmically diverging negative 
entropy. 
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0(0) 
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2(1) 
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t 0(0) 1 
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Figure 2: An extended non- vortex defect that corresponds to a well-defined set of Potts 
states (first label) but for which the heights mod 6 (second label, in parentheses) cannot 
be lifted uniquely everywhere from Z/6Z to Z, since the height at the center could corre- 
spond to either +3 or —3. Equivalently, there are two flux-conserving ways of assigning 
triple arrows to the unsatisfied bonds (encircling the center plaquette clockwise or anti- 
clockwise) . Nevertheless, sufficiently far away from the defect there is no such ambiguity 
in lifting to Z. 
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Figure 3: £(/?)/e 2/3 versus (3. Error bars are one standard deviation, and are most likely 
overestimates |10|, Section 4.1.1]. 
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Figure 4: logf (/?) - 2/3 versus /3e~ f ' '. Straight line is logf (/?) - 2/3 = -2.11 - 6.65/3e" 
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Figure 5: logf (/?) -2/3+6.65/3e"0 versus /3e _/J . Straight line is logf -2/5+6.65/?e 
-2.11. 
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Figure 6: log£( / 3)-2/3+6.65 / 3e- /3 versus (3e' 2(B . Straight line is log£(/3)-2/?+6.65/3e 
-2.12 + 123(3e~ 213 . 
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Figure 7: logx(/?)-f P versus /Je"". Straight line is logx(P)-fP = -2.52-ll.82/5e"' 3 . 
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Figure 8: logx(/3) - fP + H.82/3e _/3 versus fie' 13 . Straight line is \ogx(P) - f P + 
11.82/5e" /3 = -2.52. 
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Figure 9: logx(P) ~ fP + H.82/3e _/J versus (3e' 2/3 . Straight line is log 
11.82/3e _/J = -2.54 + 232 / 3e~ 2/3 . 
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